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Abstract 

We use polynomial truncations of the Fourier transform of the local mea- 
sure to calculate the connected g-point functions of Dyson's hierarchical model 
in the broken symmetry phase. We show that accurate values of the connected 
1, 2 and 3 point functions can be obtained at large volume and in a limited 
range of constant external field coupled linearly to the field variable. We in- 
troduce a new method to obtain the correct infinite volume and zero external 
field extrapolations. We extract the leading critical exponents and show that 
they obey the scaling and hyperscaling relations with an accuracy ranging 
from 10~ 5 to 5 x 10~ 3 . We briefly discuss how to improve the method of 
calculation. 
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I. INTRODUCTION 



Spontaneous symmetry breaking plays a fundamental role in our understanding of the 
mass generation mechanism of the elementary particles. One of the simplest field theory 
model where it is observed is scalar theory Despite its simplicity, there exists no known 
analytical method which allows one to elucidate quantitatively all the dynamical questions 
which can be asked about scalar field theory in various dimensions. From a sample of the 
recent literature on scalar field theory, one can see that the Monte Carlo method is a popular 
tool to settle questions such as the existence of non-perturbative states M, large rescaling 
of the scalar condensate or Goldstone mode effects |§. 

The Monte Carlo method allows us to approach quantum field theory problems for which 
there are no known reliable series expansions. The main limitations of the method are the 
size of the lattice which can be reached and the fact that the errors usually decrease like t -1 / 2 , 
where t is the CPU time used for the calculation. If, in the next decades, a better knowledge 
of the fundamental laws of physics has to rely more and more on precision tests, one should 
complement Monte Carlo methods with new computational tools which emphasize numerical 
accuracy. 

This motivated us to use "hierarchical approximations" as a starting point, since they 
allow a more easy use of the renormalization group (RG) transformation. Examples of hier- 
archical approximations are Wilson's approximate recursion formula Q or the hierarchical 
model ||. In the symmetric phase, we have found || that polynomial truncations of the 
Fourier transform of the local measure provide spectacular numerical accuracy, namely, var- 
ious types of errors decrease like e _j4t ", for some positive constant A of order 1 when t is 
measured in minutes of CPU time and 0.5 < u < 1. In particular, t only grows as the loga- 
rithm of the number of sites L D and the finite-size effects decay like L~ 2 when L (the linear 
size) becomes larger than the correlation length. This method of polynomial truncations 
was used to calculate the critical exponent 7 in the symmetric phase for the hierarchical 
model with estimated errors of the order of 10~ 12 . The result was confirmed by calculating 
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the largest eigenvalue of the linearized RG about the accurately determined non-trivial fixed 
point p. 

Thanks to the polynomial approximation, very accurate information can be encoded in 
a very small set of numbers. In the symmetric phase, this approximation is numerically 
stable when the number of sites becomes arbitrarily large and the high-temperature fixed 
point is reached. On the other hand, in the broken symmetry phase, numerical instabilities 
appear after a certain number of iterations following the bifurcation, and it is not possible 
to completely get rid of the finite size effects with the straightforward procedure used in the 
symmetric phase. This issue was briefly discussed in section III.E of Ref. J7). 

In this paper, we analyze the numerical instabilities of the low-temperature phase in a 
quantitative way. We show that in spite of these numerical instabilities, it is possible to 
take advantage of the iterations for which the low-temperature scaling is observed to obtain 
reliable extrapolations of the magnetization, first to infinite volume at non-zero external field 
and then to zero external field. We then present a more pratical method of extrapolation 
which we apply to calculate the connected g-point functions at zero momentum G^(0) for 
q = 1, 2 and 3. Finally, we use these calculations to extract the leading critical exponents 
and we check the hyperscaling relations among these exponents. 

The paper is organized as follows. In section II we show how to construct recursively 
the generating function for the G°(0) when a magnetic field is introduced. In section III, we 
review the scaling and hyperscaling relations among the critical exponents and explain how 
they should be understood in the case of the hierarchical model. Hyperscaling || usually 
refers to scaling relations involving the dimension explicitly. Dyson's hierarchical model 
has no intrinsic dimensionality but rather a continuous free parameter usually denoted by 
c introduced in section II, which controls the decay of the interactions among blocks of 
increasing sizes. This parameter can be tuned in order to insure that a massless field 
has scaling properties that can be compared with those of nearest neighbor models in D- 
dimensions. In the past we have chosen the parametrization c = 2 1_2 / D , however this is not 
the only possible one. In section III. C, we show that a more general parametrization of c, 
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(which includes 77) combined with linear arguments yields predictions that are identical to 
the conventional predictions obtained from scaling and hyperscaling. We want to emphasize 
that the main prediction of the linear theory - that can be interpreted as a hyperscaling 
relation - can be expressed in terms of c only and is given in general by Eq. ( |3.15| ). For 
c = 2 1 / 3 , this general equation together with the accurate result of Ref. [0] implies 

7 3 = 1.29914073... x (5g-6)/4, (1.1) 

where 7 9 is the leading exponent corresponding to the connected g-point function. 

We then proceed to verify the predictions of Eq. (|1.1|) by doing actual calculations at 
various values of the inverse temperature (3 near criticality. This is a rather challenging task 
because as one moves away from the unstable fixed point, in the low-temperature side, rapid 
oscillations appear in the Fourier transform of the local measure and the polynomial approx- 
imation ultimately breaks down. This is the cause of the numerical instabilities mentioned 
above. As a consequence, a relatively small number of iterations can be performed with a 
reasonable accuracy in the low-temperature phase. This is explained in section IV where we 
also show that the number of numerically accurate iterations in the low-temperature phase 
scales like the logarithm of the degree of the polynomial. For the calculations discussed 
later in the paper, we have used a polynomial truncation of order 200. With this choice, the 
number of iterations where an approximate low-temperature scaling is observed is slightly 
larger than 10. Since for Dyson's hierarchical model the number of sites is halved after 
each iteration, it means roughly speaking that in correlation length units we can only reach 
volumes which are 2 10 ~ 10 3 . If we use the D = 3 interpretation of c = 2 1//3 , this means 
that the linear size, denoted by L, which can be reached safely are at most 10 times the 
correlation lengths. 

Despite this limitation, the magnetization reaches its infinite volume limit with clearly 
identifiable L~ 2 corrections provided that the external magnetic field is not too large (oth- 
erwise the polynomial approximation breaks down) or not too small (otherwise a linear 
analysis applies and there is no spontaneous magnetization). The exact intermediate range 
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of the magnetic field for which the connected g-point functions reach an infinite volume limit 
with the characteristic L~ 2 corrections is discussed in section V. In this intermediate range, 
two methods of extrapolation can be used. The first is the standard one which consists in 
extrapolating to infinite volume at fixed external field and then to zero external field. On 
the other hand, within the intermediate range of magnetic field mentioned above, the mag- 
netization at finite volume can be fitted very accurately with a straight line which provides 
an extrapolation to zero magnetic field. This extrapolation has no physical meaning but it 
also reaches an infinite volume limit with L~ 2 corrections when the volume increases. This 
limit coincides with an accuracy of 6 significant figure with the limit obtained with the first 
method; in other words within the estimated errors of the calculation. The second procedure 
is much more practical because it does not require any overlap among the acceptable regions 
of magnetic field when the volume increases. The second method will be used to calculate 
the higher point functions. 

Proceeding this way, we calculate the connected g-point functions at zero momentum 
Gq(0), for q = 1, 2 and 3 and for various values of the inverse temperature (3. The results 
are reported in section VI. The critical exponents are then estimated by using a method 
discussed in Ref. [[7J where we selected a region of f3 for which the combined effects of the 
errors due to subleading corrections and the numerical round-off could be minimized. Using 
linear fits within this limited range of f3, we found exponents in agreement with the prediction 
of hyperscaling given in Eq. (|3.15|) with an accuracy of 10 -5 for the magnetization, 4 x 10~ 5 
for the susceptibility and 5 x 10~ 3 for the 3-point function. As far as the first two results are 
concerned, the accuracy compares well with the accuracy that can usually be reached with 
a series analysis or the Monte Carlo method. Nevertheless, there is room for improvement: 
one should be able to "factor out" the rapid oscillations in the Fourier transform of the local 
measure and treat them exactly. This is discussed briefly in the conclusions. 
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II. INTRODUCTION OF A MAGNETIC FIELD 



In this section, we describe Dyson's Hierarchical Model |Jl(J coupled to a constant 
magnetic field. All calculations are performed at large but finite volume. The total number 
of sites denoted 2 nmax . We label the sites with n max indices x nmax , ...,xi, each index being 
or 1. In order to visualize this notation, one can divide the 2 nmax sites into two blocks, each 
containing 2 rimax ~ 1 sites. If x nmax = 0, the site is in the first box, if x nmax = 1, the site is 
in the second box. Repeating this procedure n max times (for the two boxes, their respective 
two sub-boxes , etc.), we obtain an unambiguous labeling for each of the sites. 

The non-local part of the action (i.e. the "kinetic term") of Dyson's Hierarchical model 
reads 

Q Umax £ 

Skin = - ^ E (4)" S ( E 0(*nm«.-.*l)) 2 ■ ( 2A ) 

71=1 "^ n max i&n-\-l %n j ■ ■ ■ ,^1 

The index n, referred to as the 'level of interaction' hereafter, corresponds to the interaction 
of the total field in blocks of size 2 n . The constant c is a free parameter which describes 
the way the non-local interactions decay with the size of the blocks. We often use the 
parametrization 

c = 2 l - 2 l D , (2.2) 

in order to approximate D-dimensional models. This question will be discussed later (see 
Eq. (gig) for a generalization of Eq. Q). 

A constant external source H, called "the magnetic field" later, is coupled to the total 
field. This can be represented by an additional term in the action 

S H = -H 0(zn m a*>-,si) ■ ( 2 - 3 ) 

x rimax v^l 

However due to the linearity of the coupling, e~ SH factorizes into local pieces and this 
interaction can be absorbed in the local measure. The field 4>( Xnmax ,...,x 1 ) is integrated over 
with a local measure 
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Wo(0,#)ocWo(0)e H *, (2.4) 

where Wq(4>) is the local measure at zero magnetic field. For simplicity, we use the convention 
that if the magnetic field does not appear explicitly in an expression (e. g. , Wo(<f>)) the 
quantity should be understood at zero magnetic field. The constant of proportionality 
refers to the fact that we require both Wq(4>, H) and Wq(4>) to be normalized as probability 
distributions. Since we are interested in universal properties, we will use a single local 
measure, namely the Ising measure, Wo(<f>) = 5(<f) 2 — 1). Numerical experiments in Ref. |7j 
show that the universal properties are very robust under changes in the local measure. 

At H = 0, the recursion relation corresponding to the integration of the fields in boxes 
of size 2, keeping the sum of the 2 fields in each box constant reads 

w n+ M = ^±V/ 2 ( C /4r+v J d ^w n {^l)w n {^l) , (2.5) 

where C n+ \ is a normalization factor which will be fixed in order to obtain a probability 
distribution. Introducing the Fourier representation as in Refs. |||7|] 

W n (<j>) = J $Le**tf n (k) , (2.6) 

and a rescaling of the "source" k by a factor (c/4) 1 / 2 at each iteration , 

R n {k) = W n {k{c/A) n / 2 ) , (2.7) 

the recursion relation becomes 

Rn+i(k) = C n+l ex P {- l -f3^){R n {k{c/A) l l 2 )) 2 . (2.8) 

We fix the normalization constant C n is such way that R n (0) = 1. R n (k) has then a direct 
probabilistic interpretation. If we call M n the total field 4>x inside blocks of side 2™ and 
(...)„ the average calculated without taking into account the interactions of level strictly 
larger than n, we can write (at H = 0) 

Rn{k) = E [ -7Ar((Mn) 2q )n(c/Ar ■ (2.9) 
g=0 \ A( D- 
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The introduction of the magnetic field is a very simple operation. The basic equation 
reads 

W n (<p, H) oc W n (<j>)e H + . (2.10) 



This can be seen in many different ways. One of them is to use Eq. ( |2.5| ) and realize that 
the 0' drops out of the magnetic interactions. Another one consists in realizing that due 
to the linearity, one can split Eq. (|2.3|) into sum over boxes of any desired size. In Fourier 
transform, this implies that 

W n {k,H) oc W n (k + iH) . (2.11) 

The normalization factor is fixed by the condition W n (0,H) = 1 which guarantees that 
W n (<fi, H) is a probability distribution and that W n (k,H) generates the average values of 
the positive powers of the total field. More explicitly, 

W n (k,H) = W ^ + m) (2.12) 
W n {iH) 

From a conceptual point of view, as well as from a practical one, it is easier to deal with 
the rescaled quantity R n {k). Near the fixed point of Eq. (|2.8|), we have the approximate 
behavior 

((M n ) 2 % oc (4/cr , (2.13) 
In terms of the rescaled function, we can rewrite Eq. ( [2.131 ) as 

R n ^H(4/c)^) ~ ^ ^r ((Mn) )nAC,A) ■ (2 ' 14) 

The connected Green's functions can be obtained by taking the logarithm of this generating 
function. 
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III. ABOUT HYPERSCALING 



A. General Expectations 

The main numerical results obtained in this paper are the calculations of the critical 
exponents corresponding to the singularity of the connected g-point functions for q = 1, 2 
and 3. For defmitiness we use the notation 

^(0)oc(/3-/3 c )-^ , (3.1) 

for the leading singularities in the low-temperature phase. We assume that the reader is 



familiar with the commonly used notations |TT| for the critical exponents. For q = 1, we have 



7i = —j3 which should not be confused with the inverse temperature. After this subsection, 
we keep using the notation (3 for the inverse temperature. For q = 2, we have 72 = 7'. If 
one assumes that the scaled magnetization Mj (T — T c )@ is a function of the scaled magnetic 
field H/ (T — T C ) A only, one obtains that 

7 9+ i - l q = A , (3.2) 

for any q. The exponent A is often called the gap exponent and should not be confused 
with the exponent associated with the subleading corrections to the scaling laws. 

In general, there exists 7 relations among the 10 critical exponents a, a' ,(3 ,7 ,7', 
A , 5 , v , v' , and r], in which the dimension of the system does not enter explicitly. These 



are the so-called scaling relations [TT[] which stimulated the development of the RG method. 
Their explicit form is 

a = a' (3.3) 

7 = 7' (3.4) 

v = v' (3.5) 

a + 2/5 + 7 = 2 (3.6) 

A = /3 + 7 (3.7) 



A = (35 (3.8) 
J=(2-v>- (3.9) 

Eq. ( |3.7| ) can be seen as an obvious version of Eq. ( |3.2| ) for q = 1 , but has also a non-trivial 
content summarizing Eq. ( |3.2[ ) for all the higher q. 

In addition there exists one relation where the dimension enters explicitly, for instance: 

Dv = 2-a. (3.10) 

Other relations may be obtained by combining Eq. ( p.lOj ) with the scaling relations. Pro- 
ceeding this way, we obtain a relation of relevance for the rest of the discussion, namely 

H 2(2-r]) 1 K ' 

The relations involving the dimension explicitly are usually called hyperscaling relations [|12| . 



A mechanism leading to a possible violation of hyperscaling (dangerous irrelevant variables) 
is explained in appendix D of Ref. ||. If the 8 relations hold, we are left with only two 
independents exponents, for instance 7 and r\. 

Combining the hyperscaling relation ( 3.11| ) and the scaling relations ( |3.4j ) and ( |3.7| ), we 



obtain 



7g = 7 +(g-2)A (3.12) 
= 7 [-2£> + q(D + 2 - »j)]/(4 - 2r]) 



B. The Hierarchical Model (HT case) 

In the case of the hierarchical model, the exponents 7 9 of the high-temperature (HT) 
phase (so for q even) can be estimated by using the linearized RG transformation. Since this 
subsection is the only part of this article where we will consider the high-temperature phase, 
we have not found useful to introduce special notations for 7 q in this phase. When f3 c — (3 is 
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small, the linearized RG transformation can be used for approximately n* iterations, with 
n* defined by the relation 

\(3 - (3 c \\ n * = 1 , (3.13) 

where A is the largest eigenvalue of the linearized RG transformation. After the transient 
behavior has died off and until n reaches the value n*, we are near the fixed point and R n (k) 
does not change appreciably. Remembering that the field is rescaled by a factor (c/4)2 at 
each iteration (see Eq. (|2.8|) ), we obtain the order of magnitude estimate for G q (0) after n* 
iterations: 

G c q {0) « 2~ n *(4/c) 9n */ 2 . (3.14) 

For n larger than n*, the non-linear effects become important. The actual value of G q (0) 
may still change by as much as 100 percent, however the order of magnitude estimate of Eq. 
Q3.14j ) remains valid. This transition has been studied in detail in Ref. Jl3| in a simplified 



version of the model. Eliminating n* in terms of (3 C — (3, we obtain the value of the leading 
exponents 

lq = 7 [(g/2)ln(4/c) - ln2]/ln(2/c) , (3.15) 

with 

7 = ln(2/c)/lnA . (3.16) 

This relationship has been successfully tested || in the symmetric phase for q = 4 and 
4/c = 2 5 / 3 . 

C. About Dimensionality 



We will now show that Eq. ( |3.15| ) is compatible with the general relation of Eq. fl3.12| ) 



provided that we relate c to a parameter D which can be interpreted as the dimension of 
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a nearest neighbor model approximated by the hierarchical model. We introduce a linear 
dimension L such that the volume L D is proportional to the total number of sites 2 n . From 

Loc2 n / D , (3.17) 

we can in general relate c and D by assuming a scaling of the total field 

(M^) n ocL (D+2 -" )9 , (3.18) 

From comparison with Eq. ( |2.13| ) this would imply that 

(4/c) = 2 {D+2 -^ )/D . (3.19) 



Substituting in Eq. (|3.15 ), we reobtain the general Eq. ( |3.12| ) 



Since in the infinite volume limit, the kinetic term is invariant under a RG transformation, 
we have chosen in the past to use Eq. ( |3.19j ) with rj = 0. This is our conventional definition 



of c given in Eq. (|2.2j ). This is the same as saying that when we are near the fixed point, 
the total field in a box containing 2 n sites, scales with the number of sites in the same way 
as a massless Gaussian field. This obviously implies that in the vicinity of a Gaussian fixed 
point the total field scales exactly like a massless Gaussian field in D dimension. On the 
other hand, an interacting massless field will also scale like a free one, which is not a bad 
approximation in D = 3. This is an unavoidable feature which will need to be corrected 
when one tries to improve the hierarchical approximation. 

We emphasize that this interpretation has no bearing on the validity of the calculations 
performed. What matters in our calculation is the value of 4/c. In the following, we have 
used 4/c = 2 5 / 3 , which can be interpreted either as D = 3 and r] = or, for instance, as 
D = 2.97 and rj = 0.02. 



D. The Low- Temperature Case 



The extension of the argument for odd and even values of q in the broken symmetry 
phase is somehow non-trivial. Since we need to take the infinite volume limit before taking 
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the limit of a zero magnetic field, we need some understanding of the non-linear behavior. 
Some aspects of the non-linear behavior are discussed in section |V|. In the following, we 
will show numerically that Eq. ( |3.15| ) holds in good approximation in the broken symmetry 



phase for 4/c = 2 5 / 3 . With this choice of 4/c and the corresponding value of 7 calculated 
in Ref. |7|], Eq. ( |3.15| ) implies Eq. ( |1 . 1| ) given in the introduction. The verification of this 



relation for q—1, 2 and 3 is the main numerical result discussed in the following chapters. 

IV. POLYNOMIAL TRUNCATIONS 

In the following we will exclusively consider the case of an Ising measure 

R (k) = cos(k) . (4.1) 

This restriction is motivated by accurate checks of universality based on calculations with 
other measures. Given that Rq can be expanded into a finite number of eigenfunctions of 
exp(— 1/3 Jp-), one can in principle obtain exact expressions for the next R n (k), for instance 

1 + gco#) 
Rl W = i + e /3c/ 2 • ( 4 - 2 ) 

One can in principle repeat this procedure. At each iteration, one obtains a superposition 
of cosines of various frequencies. For a given numerical value of c, n iterations of this exact 
procedure requires to store 2 ra_1 + 1 numerical coefficients. The memory size thus scales 
like 2 n , while the CPU time scales like 4 n . If (3 differs from f3 c by 10~ 10 , one needs at least 
80 iterations in order to eliminate the finite-size effects. Such a calculation using the exact 
method can be ruled out by practical considerations. 

We will thus try to extend the approximate methods that we have used successfully in 
the symmetric phase J7|, where the function R n {k) was calculated using finite dimensional 
approximations || of degree I max'- 

R n (k) = 1 + a nA k 2 + a nt2 k 4 + ... + a nUax k 2lma * . (4.3) 
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After each iteration, non-zero coefficients of higher order {dn+i,i max +\ e ^ c - ) are obtained, but 
not taken into account as a part of the approximation in the next iteration. The recursion 
formula for the a n ^ m reads |J : 

E|—(E P+g= z«n, P an, g )[(2/)!/(/ - m)\(2m)\](c/AY[-(l/2)P] l - m 
ttn+1 ' m ~ E!=r(E p+9 =za n , p a ni9 )[(20!/n]( C /4)'[-(l/2)/3]' ^ 1 ' ' 

The method to identify f3 c has been discussed in detail in Ref. |J and consists in finding 
the bifurcation in the ratio a n+ \^/ a n ^. In the following, we simply call this quantity "the 
ratio". If j3 < /3 C , the ratio drops to c/2 for n large enough. In this case, the numerical 
stability of the infinite volume limit is excellent and allows extremely accurate determination 
of the renormalized quantities. If (3 > @ c , the ratio "jumps" suddenly a few iterations 
after n* is reached and stabilizes near the value c, corresponding to the low-temperature 
scaling. This is seen from Eq. (|2.13| ). Since (M 2 ) n grows like 4 n , as one expects in the 
low-temperature phase, and remembering that there is a rescaling of c/4 at each iteration, 
the coefficient of k 2 grows like c n . This implies a ratio equal to c. In our calculation, 
c = 1.25992 . . .. Unfortunately, the number of iterations where the low-temperature scaling 
is observed is rather small. Subsequently, the ratio drops back to 1. As we shall explain at 
length, this is an effect of the polynomial truncation. The length of the "shoulder" were the 
low-temperature scaling is observed increases if we increase l ma x- This situation is illustrated 
in Figure 0. 
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FIG. 1. The low-temperature shoulder at = (3 C + 10 1 for l max = 200 (empty circles) and 
l-max = 80 (filled circles) as a function of l max . 

No matter how large l max is, for n large enough, the ratio eventually drops back to 1. 
This reflects the existence of a stable fixed point for the truncated recursion formula. The 
values a* of ai at this fixed point for various l max are shown in Figure 0. 
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FIG. 2. Value of for the "false" low-temperature fixed points for the Ising case in 3 dimensions 
for I = 1 (circles), I = 2 (filled stars), I = 3 (empty stars) 



We see a clear evidence for a dependence of the form 
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a* cc (l max y . (4.5) 

This means that the stable fixed point is an effect of the polynomial truncations and has no 
counterpart in the original model. 

It is possible to evaluate the value of n for which the low-temperature shoulder ends. A 
detailed study shows that for n large enough, we have in good approximation 

R n {k) ~ cos(A^c n/2 A;) , (4.6) 

where M. is the magnetization density in the infinite volume limit. If we assume that 
R n (k) is exactly as in Eq. Q4.6D , then we can use the basic recursion formula (|2.8| ) in 
order to obtain the corresponding R n+ i(k). Using 2 x cos 2 (x) = 1 + cos(2a;), we can re- 
express (i? n (fc(c/4) 1 / 2 )) 2 as a superposition of eigenfunctions of the one- dimensional Lapla- 
cian. When the exponential of the Laplacian in Eq. (|2.8| ) acts on the non-constant modes 



it becomes exp(/3A / t 2 c n+1 /2). In the polynomial truncation of the recursion relation, this 
exponential is replaced by l max terms of its Taylor expansion. This approximation is valid 
if the argument of the exponential is much smaller than l m ax- Consequently, we obtain that 
the polynomial truncation certainly breaks down if n is larger than n b such that 

n b + 1 ~ [ln(2//3) - ln(.M 2 ) + ln(Z mM )]/mc . (4.7) 

If the estimate of Eq. Q3.14|) extends to the low-temperature phase, one realizes that the 
second term of (f4.7|) is roughly n* while the third term stands for the length of the peak and 
the shoulder. Plugging the approximate values 1.1 for (3 and 0.7 for M. (see section |V|), we 
obtain rif, = 23 for l max = 80 and n b = 27 for l max = 200. A quick glance at Figure [j], shows 
that these estimates coincide with the first drastic drops of the low-temperature shoulder. 

One can in principle extend indefinitely the low-temperature shoulder by increasing l ma x- 
However, the CPU time t necessary for n iterations of a quadratic map in dimension l max 
grows like 

t oc n(l max ) 2 . (4.8) 
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As we will show in section [V|, the finite-size effects on Gg(0) are of the order (c/2)" s where 
n s is the number of points on the shoulder. This behavior has been demonstrated in 
the high-temperature phase and we will see later that it also applies in the low-temperature 
phase. From the previous discussion n s « lnZ max /lnc. This implies that the finite-size effects 
£ are of the order 

8<x{l max )^l l " c . (4.9) 

Using Eq. ( |4.8| ) and the value of c expressed in terms of D according to Eq. ( |3.19| ) with 
i] = 0, we obtain 

Ext- 1 " -** . (4.10) 

In particular, for the value 4/c = 2 5//3 used hereafter, the errors decrease like t" 1 . Con- 
sequently, we should try to modify the method in such a way that the rapidly oscillating 
part of R n (k) is treated without polynomial approximations. This possibility is presently 
under investigation. One can nevertheless obtain results with an accuracy competing with 
existing methods by using the finite data on the short shoulder in order to extrapolate to 
the infinite volume limit result. This procedure is made possible by the rather regular way 
the renormalized quantities approach this limit. 



V. THE EXTRAPOLATION TO INFINITE VOLUME 



A. Preliminary Remarks 

There is no spontaneous magnetization at finite volume. This well-known statement can 
be understood directly from Eq. (|2.14 ). As explained at the beginning of section at 



finite n, R n (k) is simply a superposition of cosines with finite positive coefficients provided 
that (3 is real. However if f3 is complex, these coefficients have singularities. This comes 
from the normalization factor, needed when we impose the condition i?„(0) = 1, which has 



zeroes in the complex plane. The behavior of these zeroes has been studied in Ref. [0] for n 
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between 6 and 12. As the volume increases, these zeroes "pinch" the critical point. However 
at finite n, there are no zeroes on the real axis. In conclusion, at real j3 and finite n, R n {k) 
is an analytical function of k. For any given n, we can always take the magnetic field ff 
small enough in order to have 

|ff(4/c) n/2 | < 1 . (5.1) 



If we express c in terms of the linear dimension using Eqs. fl2-2|) and (|3.17|) this translates 
into 

|ff | < L-( D+2 V 2 . (5.2) 

Given the analyticity of R n (k), one can then use Eq. (|2.14j) in the linear approximation. 
In this limit, 

(M n ) n ~ -2a n)1 ff(4/c) n , (5.3) 

and the magnetization vanishes linearly with the magnetic field. 

On the other hand, for any non-zero ff , no matter how small its absolute value is, one 
can always find a n large enough to have |ff(4/c) n / 2 | > 1. The non-linear effects are then 
important and Eq. Q5.3p does not apply. In addition it is assumed (and will be verified 
explicitly later) that when such a n are reached, the value of the G^(0) stabilizes at an 
exponential rate. One can then, first extrapolate at infinite volume for a given magnetic 
field, and then reduce the magnetic field in order to extrapolate a sequence of infinite volume 
limits with decreasing magnetic field, towards zero magnetic field. Again, this procedure 
requires some knowledge about the way the second limit in reached. In the case considered 
here (one scalar component), the limit is reached by linear extrapolation. In systems with 
more components, the Nambu-Goldstone modes create a square root behavior [[15] 

M(T < T c , ff > 0) = M(T, + ) + CH 1 / 2 . (5.4) 

which has been observed for 0(4) models using Monte Carlo simulations [Q. We now discuss 
the application of the procedure outlined above in the simplest case. 



B. Calculation of the magnetization 



In this subsection we discuss the calculation of the infinite volume limit of the magneti- 
zation. The magnetization density at finite volume is defined as 

M n (H) = . (5.5) 

We call it "the magnetization" when no confusion is possible. For definiteness, we have 
chosen a special value (3 = fic + lO" 1 and calculated the magnetization by plugging numerical 
values of H in Eq. Q2.14 ) and expanding to first order in k. The results are shown in Fig. 



H for n = 17 and l max = 200. 
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FIG. 3. log 10 (M) versus log 10 (#) at n = 17 for (3 = C + 10 -1 . 

As one can see, we have three different regions. The first one (I) is the region where 
the linear approximation described above applies. For the example considered here, the 
linearization condition \H{A/ c) n ^ 2 \ 1 translates into log 10 (i7) <C —4.3. This is consistent 
with the fact that the linear behavior is observed below -5. The third part (III), is the region 
where the polynomial approximation breaks down. Given the approximate form given in Eq. 
Q16D , this should certainly happen when \H(A/c) n ^ 2 \ ~ l ma x- This means \og w (H) w —2.0 
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in our example. On the figure, one sees that for log l0 (H) rs —2.4, the magnetization drops 
suddenly instead of reaching its asymptotic value at large H, namely A4 — 1. Finally, the 
intermediate region (II) is the one which contains the information we are interested in. 

As advertized, we will first take the infinite volume limit of the magnetization at non- 
zero magnetic field and then extrapolate to zero magnetic field. We need to understand how 
the second region shown in Fig. ^| changes with n. From the above discussion, region II is 
roughly given by the range of magnetic field 

- (7i/2)log 10 (4/c) < log 10 (#) < \og 10 (l max ) - (n/2)log 10 (4/c) . (5.6) 

In the log scale of Fig. |3|, the width of region II is at most log 10 (l max ) which is approximately 
2.3 in our sample calculation. Region II shifts by — (l/2)log 10 (4/c), approximately 0.25 in 
our sample calculation, at each iteration. In addition, the whole graph moves slightly up at 
each iteration in a way which is better seen using a linear scale as in Fig. |J 
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200 and (3 = C + 10 
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As one can see, the regions II of seven successive iterations do not overlap. Consequently 
1.5 is a more realistic estimate than the previously quoted bound 2.3 for the average width 
of region II. 
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The fewer iterations we use to extrapolate to infinite volume, the broader the range of the 
magnetic field can be. We have compared 5 sets of 4 iterations well on the low-temperature 
shoulder starting from the set (17, 18, 19, 20) up to the set (21, 22, 23, 24). From our 
experience in the symmetric phase we have assumed that the finite-size effects could be 
parametrized as 

M n = Moo-AxB n . (5.7) 

This relation implies that 

\og w (M n+1 - M n ) = A + nx \og w (B) , (5.8) 

where A = \og lQ (A) + log 10 (l — B). The value A and log 10 (S) can be obtained from linear 
fits. For four successive iterations, this will give us 3-point fits for the infinite volume limit 
at fixed value of H ^ 0. In all fits performed, we found B ~ 0.63, which is compatible with 
the (c/2) n decay of the finite size effects found in the symmetric phase ||. In terms of the 
linear dimension L introduced in Eq. ( |3.17| ), this corresponds to finite-size effects decaying 



like L~ 2 . If the parametrization of Eq. ( |5.7| ) was exact, the value of Ai n + A x B n would 
be independent of n and equal to M.^. In practice, variations slightly smaller than 10~ 6 
are observed. We have thus taken an average over these values in order to estimate .Moo at 
fixed H . The results for the first set are shown in Fig. |5] for various values of H . The linear 
behavior allows an easy extrapolation to H = 0. 
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Magnetic Field H 



FIG. 5. Moo extrapolated from the data for n= 17, 18, 19 and 20, versus the magnetic field for 
/ ma:r = 200 and/3 = /3 c + 10- 1 . 

We have repeated this procedure for the four other set of successive values defined pre- 
viously and obtained the H = extrapolations: 



Set M»=° 

1 0.7105296 

2 0.7105349 

3 0.7105376 

4 0.7105380 

5 0.7105382 



Averaging over these five values, we obtain 



MH =0 = 0.710536 =f 3 x 1(T 6 . (5.9) 

It may be argued that the values coming from sets involving larger values of n are better 
estimates because the finite size effects are smaller for those sets. 
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We have repeated this type of calculation with sets of 5 successive iterations and a 
correspondingly narrower range of magnetic field and found results compatible with the 
estimate given by Eq. (|5.9|) . 



C. The Susceptibility 

We now consider the calculation of the connected susceptibility (two-point function). By 
using the previous notation, we can express it as 

(M n )l H - (MX,h 



Xn(H) 2 „ 

(6? - 2 x b 2 ) 

1 " (5.10) 



2 n 

where 

R n (k + iH(A/c)^) " 

The dependence on H of the h n is implicit. 

In order to extrapolate the susceptibility to infinite n, one has to determine the range 
of the magnetic field for which the scaling (M n )^ H — (M^ H ) n cx 2 n holds. When this is the 
case, the ratio Xn+i/Xn — 1- The range of values of the magnetic field for which this scaling 
is observed is analogous to "region II" introduced in the previous subsection, and we will 
use the same terminology here. The ratios of the susceptibility at successive n are shown 
for various values of H in Fig. 13. 
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FIG. 6. The successive ratios of the susceptibility for (3 = (3 C + 10 . 

One observes that the range where the desired scaling is observed shrinks when n in- 
creases. For each successive iteration, the ratio of susceptibility has a "upside-down U"- 
shape. The values of H for which the ratio starts dropping on the left are equally spaced 
and can be determined by linearization as before. On the other side of the upside-down 
U, dropping values of the ratio signal the breakdown of the polynomial truncation. This 
occurs at smaller values of H than for the magnetization, making the region II smaller. A 
theoretical estimate of the lower value of H for which this happens requires a more refined 
parametrization than the one given in Eq. (|4.6|) . In order to get a controllable extrapolation, 
we need at least 4 successive values of \ n (to get at least 3-point fits for the logarithm of the 
differences). This is unfortunately impossible: the region II of three successive upside-down 
U have no overlap as one can see from Fig. |6|. Similar results are obtained by plotting the 
susceptibility versus the magnetic field as shown in Fig. [7| 
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FIG. 7. The susceptibility \ n versus the magnetic field H at different n for (3 = /3 C + 10 . 

One sees that the regions II (where the susceptibility can approximately be fitted by 
a line which extrapolates to a non-zero value at zero H) do not overlap for 4 consecutive 
iterations. 

D. An Alternate Method 



In the previous discussion, we have observed a linear behavior for the region II of the 
magnetization and the susceptibility. This linear behavior can be used to obtain extrapo- 
lations to non-zero values of these quantities at zero magnetic field. These values have no 
physical interpretation. We denote them by A4^°, the arrow indicating that the quantity 
is a mathematical extrapolation and not "the spontaneous magnetization at finite volume" . 
They reach an asymptotic value at an exponential suppressed rate when n increases, just as 
in Eq. (|5.7|). This is illustrated in Fig. |8L 
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FIG. 8. log 10 (7W^° - M^°) versus n for l max = 200 and (3 = f3 c + 10" 1 . 

Using a linear fit to fix the unknown parameters A and B in Eq. |5]7], and averaging the 
■M-H^ + AB n over n, we obtain 



M" =0 = 0.710537 =f 2 x 10~ 6 



(5.12) 



which is consistent with the result obtained with the standard method. 

Roughly speaking, the lines of region II move parallel to each other when n increases 
and it is approximately equivalent to first extrapolate to zero H, using the linear behavior 
in region II, and then to infinite n rather than the contrary. If the two limits coincide, the 
second method has a definite practical advantage: all we need is a small part of region II for 
each n, no matter if it overlaps or not with the region II for other n. So in general, it allows 
to use more iterations to get better quality extrapolations. The fact that the values of the 
magnetization obtained with the two methods coincide with 5 significant digits is a strong 
indication that the two procedure are equivalent. For the susceptibility and higher point 
functions, we do not have an independent check, since the alternate method is the only one 
available. However, we were able to make consistency checks such as the fact that the slope 
of the straight line used for the zero magnetic field extrapolation of the q— point function 
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coincides with the q + 1-point function. 

We can repeat the same steps for the 3-point function. The 3-point function is given by 

M 3 - 3MiM 2 + 2M 3 



G 



3 2 r 



Q h _ Q blb2 _ 2b\ 

^ > ( 5 - 13 ) 



where the dependence on H is implicit. As shown in Ref. |T7|], G3 < for H < 0. Due to the 
additional subtraction, the range where the proper low-temperature scaling is observed is 
smaller than for the suceptibility. It is not possible to repeat the same steps for the 4-point 
function which is given by, 

G\ = M 4 - 3M 2 2 - 4M!M 3 + 12M 2 - 6M? , (5.14) 

and involves one more subtraction. 

VI. ESTIMATION OF THE EXPONENTS 

We have used the method described in the previous section to calculate the value of the 
connected g-point functions at value of f3 approaching f3 c from above with equal spacings on 
a logarithmic scale. For reference, the numerical values are given in the table below. 



09 - Pc) 


Gf(0) 


GI(0) 


ci(o) 


1 


0.710537 


5.1449 


-452 


2 


0.372929 


147.75 


-4.83 x 10 6 


3 


0.181173 


3270.0 


-4.42 x 10 s 


4 


8.64639 x 10~ 2 


67534 


-3.82 x 10 11 


5 


4.10479 x 10~ 2 


1.3628 x 10 6 


-3.23 x 10 14 


6 


1.94518 x 10~ 2 


2.7276 x 10 7 


-2.72 x 10 17 


7 


9.21183 x 10~ 3 


5.4411 x 10 8 


-2.28 x 10 20 


8 


4.36138 x 10~ 3 


1.0842 x 10 10 


-1.91 x 10 23 
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9 
10 
11 
12 



2.06473 x 10~ 3 
9.77434 x 10~ 4 
4.62716 x 10~ 4 
2.19084 x 10~ 4 



2.1596 x 10 11 
4.3010 x 10 12 
8.5641 x 10 13 
1.7042 x 10 15 



-1.60 x 10 26 
-1.34 x 10 29 
-1.11 x 10 32 
-9.40 x 10 34 



The estimated errors on the values quoted above are of order 1 in the last digit for the 
first lines of the table and slowly increase when one moves down the table. For the last lines, 
the effects of the round-off errors become sizable. Otherwise, the errors are mainly due to 
the extrapolation procedure. We have checked that the numerical values of the quantities 
at finite H had reaches their asymptotic values (well within the accuracy of the final result) 
as a function of l max . 

The results are displayed in Fig. in a log- log plot. 
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FIG. 9. log 10 (G«(0)) versus log 10 (/3 - (3 C ) for q=l, 2 and 3. 

The departure from the linear behavior is not visible on this figure. In the symmetric 
phase, we know |7] that the relative strength of the subleading corrections is approximately 
-0.57(/9 c — /3) 0,43 . It is likely that a similar behavior should be present in the low-temperature 
phase. Consequently, taking into account the data on the left part of Fig. |9] will distort the 
value of the exponents. On the other hand getting too close to criticality generates large 
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numerical errors. Using a linear fit of the data starting with the 5-th point and ending with 
the 10-th point, we obtain the value of the exponents 

71 = -0.3247 

7 2 = 1.2997 (6.1) 

73 = 2.9237 . 

This can be compared with the predictions from scaling and hyperscaling given by Eq. ( |1 . 1\) 
and amount numerically to: 

7! = -0.324785 

7 2 = 1.2991407 (6.2) 
7s = 2.923066. 

Better estimates can be obtained by using the method developed in Ref. 0, where it was 
found that the combined effects of the two types of errors are minimized for 10~ 10 < \/3 c — f3\ < 
10~ 9 . This allowed estimates of 7 (in the symmetric phase) with errors of order 3 x 10~ 5 
(compared to more accurate estimates). Using 10 values between 10 -9 and 10~ 10 with equal 
spacing on a logarithmic scale, we obtain here 

71 = -0.324775 ± 2 x 10~ 5 

72 = 1.29918 ± 10" 4 (6.3) 

73 = 2.928 ± 10~ 2 . 

The errors due to the subleading corrections and the round-off errors are approximately of 
the same order in this region of temperature 0. The errors due to the subleading corrections 
are larger for larger values of \/3 c — /3| while the numerical errors are larger for smaller values 
of \ fi c — /3\ ■ We have estimated the errors due to the subleading corrections by performing the 
same calculation between 10 -8 and 10 -9 . The errors bars quoted above reflect the differences 
with the exponents obtained in this second region. 
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VII. CONCLUSIONS 



One sees clearly that our best estimates of the critical exponents (Eq. ( p.3|) ) are fully 
compatible with the predictions of hyperscaling (Eq. ( |6.2| )). The differences between the 
predicted and calculated values are 10" 5 for 71, 4 x 10~ 5 for 72 and 5 x 10~ 3 for 73. They 
fall well within the estimated errors. Since hyperscaling is a reasonable expectation, this 
also shows that the non-standard extrapolation method that we have used is reliable. As 
far as 71 and 72 are concerned, the errors bars are smaller than what can usually be reached 
using a series analysis or Monte Carlo simulation. Our result for 7! is also compatible with 
the result 0.325 obtained in Ref. |16| for the hierarchical model (for a/d = 2/3 with their 
notations) using the integral formula. 

One could in principle improve the accuracy of these calculations by increasing the size 
of the polynomial truncation. However, the efficiency of this procedure (errors decreasing 
like the inverse of the CPU time) is not compatible with our long term objectives (errors 
decreasing exponentially). The main obstruction to keep using the polynomial truncation is 
that the generating function R n (k) starts oscillating rapidly in the low-temperature phase 
making the approximation of the exponential of the Laplacian by a sum inaccurate. It is thus 
important to obtain an approximate parametrization of R n (k) in terms of eigenfunctions of 
the Laplacian. A step in this direction is made by the parametrization of Eq. ( ^4.6| ). This 
approximate analytical form needs to be improved in order to include the connected 2-point 
and higher point functions in terms of an exponential function. This possibility is presently 
under investigation. 

This research was supported in part by the Department of Energy under Contract No. 
FG02-91ER40664. 
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I. SPONTANEOUS MAGNETIZATION 



In this section, we will explain how to calculate the magnetization. We will do this in two 
different ways. One of them is to reach to infinite volume limit at different magnetic fields 
and then extrapolate to zero magnetic field and find the spontaneous magnetization. The 
second method is to obtain the spontaneous magnetization first and then take the infinite 
volume limit. 

The major difficulty which force us to go to the second method is the limitations of the 
polynomial truncation. In the low temperature phase, we have a very limited shoulder. One 
way of increasing the number of the points on this shoulder is to increase l max but this will 
not be possible in the double precision Fortran if we want to go beyond l max = 80. For this 
reason, we will use Mathematica where we can set l max = 200. This will be very important 
when we want to take the infinite volume limit first. This shoulder is shown in the Fig. ?? 
below for l max = 80 and l max = 200. 
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FIG. 1. Shoulder at (5 = /3 C + 10 -1 . l max =200 

(empty circles) and l max — 80. (filled circles) 



Before we go into the details, we need to discuss the change of the magnetization with 



1 



respect to magnetic field. We will obtain the infinite volume limit at different magnetic field 
values. The shape of the magnetization is shown in the Fig.??. 
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As one can see, we have three different region. The first one goes linearly to zero where 
the system is not magnetized and we have weak field. The second part is which we are going 
to use to reach infinite volume limit and the third one is the part where the polynomial 
truncation stops being adequate. 



II. C M ax = 200 

A. Prom infinite volume to spontaneous magnetization 

In this subsection, we will first obtain the infinite volume limit of the magnetization and 
then make extrapolation to the zero magnetic field to obtain the spontaneous magnetization. 
The region we use is the second region shown in Fig.?? above. Since we are going to take 
the limit first, we should have this region for all shoulder points. Unfortunately, it will not 
be the case. This region will shift and give us a different part for each shoulder point. In 
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this case, It will not be possible to reach to the infinite volume limit. However, there will be 
an overlapping area for 4 successive point and the fifth point will be out of this range. This 
is shown in Fig.?? below. 
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FIG. 3. Magnetization vesus magnetic field for n = 16, 17, 18, 19 for l Ti 
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200 and 



For l max = 200, our shoulder points goes from 16 to 25 and since we can only have 4 points 
in a linear region to take the limit we make sets of four points starting from 16, 17, 18, 19 
the first set, 17, 18, 19, 20 the second set and so on until 22, 23, 24, 25. On the other hand 
the points 24 and 25 are at the end of the shoulder and for these points we will not be able 
to see the second region. Then, we will work with the first 5 sets. These four point will give 
us 3-point fits for the infinite volume limit at H ^ 0. 



B. Accuracy of the 3-point fits 

We would like to know the accuracy of out three point fits. We will assume the form of 

M n = Moo ~ A x B n (1) 
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and from this relation , 

\og w (M n+1 - M n ) = A + n x log 10 (£) (2) 

where A = \og w (A) + log 10 (l — B). We listed the errors on the three-point fits for the first 
set at different magnetic fields below. 



Magnetic Field 


A 


SA 


logic (£) 
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Vfit 


0.0001221271676601683 
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0.00065146 
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0.0001346497876601683 





355549 





0110646 


-0.205466 


0.00065011 


0.000919394 


0.0001471724076601683 





354946 





0110421 


-0.205457 


0.000648786 


0.000917522 


0.0001596950276601682 





354345 





0110197 


-0.205448 


0.000647472 


0.000915663 


0.0001722176476601682 





353746 





0109975 


-0.205439 


0.000646165 


0.000913815 


0.0001847402676601683 





353148 





0109754 


-0.20543 


0.000644866 


0.000911978 


0.0001972628876601682 





352552 





0109534 


-0.205421 


0.000643574 


0.000910151 


0.0002097855076601684 





351958 





0109315 


-0.205413 


0.000642289 


0.000908334 


0.0002223081276601682 





351365 





0109098 
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0.0002348307476601682 





350773 





0108881 


-0.205395 


0.000639741 


0.000904731 


0.0002473533676601682 





350184 





0108666 


-0.205386 


0.000638478 


0.000902945 


0.0002598759876601683 





349595 





0108453 


-0.205378 


0.000637222 


0.000901168 


0.0002723986076601683 





349009 


0.010824 


-0.205369 


0.000635972 


0.0008994 


0.0002849212276601683 





348424 


0.0108024 


-0.205361 


0.000634705 


0.000897609 


0.0002974438476601681 





347854 


0.0107735 


-0.205353 


0.000633007 


0.000895207 


0.0003099664676601682 





347461 


0.0106387 


-0.205356 


0.000625088 


0.000884008 



From these values, we then calculate the .Moo at each different magnetic field value. 
Finally, we fit all these .Moo values to obtain the spontaneous magnetization in the infinite 
volume limit. This is show in the Fig.??, below. 
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Magnetic Field H 



FIG. 4. Moo versus magnetic field for n = 16, 17, 18, 19 for l max = 200 and (3 = f3 c + 10 1 and 
the extrapolation to the zero magnetic field. 



Set M%T° SM%T l 

1 0.7105296 4.6E-7 

2 0.7105349 8.2E-7 

3 0.7105376 8.5E-7 

4 0.7105380 1.7E-7 

5 0.7105382 4.5E-7 



From these data points we conclude that 

M* =0 = 0.710536 T 3.25 x 10" 6 (3) 
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C. Prom extrapolation to the infinite volume limit 

In this subsection, we would like to extrapolate to the spontaneous magnetization and 
then take the infinite volume limit and then compare our result with the one we obtained by 
using the first method above. This comparison will be handy when we use double precision 
arithmetic with l max = 80. This procedure is quite simple. We use all the possible shoulder 
points we have and then from the linear region we make an extrapolation to the zero magentic 



field. 


The results for each n are given in the table below. 


For each fit 


n 


M H =o 


SM H =o 


16 


0.7073963579884217 


4.91E-7 


17 


0.708571938893955 


1.46E-7 


18 


0.7093049288248253 


4.71E-8 


19 


0.7097633112129702 


4.38E-9 


20 


0.7100511622588243 


1.26E-9 


21 


0.7102319520698038 


1.24E-9 


22 


0.7103455954378881 


1.08E-9 


23 


0.7104170824449012 


1.08E-9 



Now, we will use these data above to estimate the infinite volume limit. We will again use 



the form of 

M n = Moo - A x B n (4) 
This will lead us a 7-point fit which is shown in the Fig.?? below. 
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/?=/? c +io 1 _ 
U= 200 - 

_i i i i i i _i 

16 18 20 22 

n 

FIG. 5. log^-M^jf] — Mn =0 ) versus n for l max = 200 and j3 = j3 c + 10 _1 and the extrapolation 
to the infinite volume limit. 

We found that a fit = 1.5 x 10~ 3 . From this linear fit of the form \og w (M„+i - M„= ) = 
A + n\og w (B) we found that log 10 B = —0.2025442151072931 =F 0.000295107922264223 
which gives B = 0.6272718307858014 ~ 0.63 and A = log 10 A + log 10 (l — B) = 
0.3087245768455599 T 0.00563802893830401. From the last relation, we find A = 
5.461756833437532. Then, all we need to do is to put these values back into the equa- 
tion above and calculate the spontaneous magnetization A4^T°. We did this for each n 
value we have the results are listed in the table below. 



, , CO 

a 1 



can 
o 



in 

CO 

I 



I 



16 0.7105341555342905 

17 0.7105401909051874 

18 0.7105395578673588 

19 0.7105377592328215 

20 0.7105369516860849 

21 0.7105366740932179 
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22 0.7105367389793958 

23 0.7105369814041257 



By using these data point, we conclude (with the simple average and standard deviation 
ralations) that 



M" =0 = 0.710537 T 1-74 x lO" 6 (5) 



D. conclusion for l max = 200 



By comparing the two results obatined by two different ways for l max = 200 and j3 = 
j5 c + 10 _1 , we immedeately see that there is 5 significant digits accuracy between them. 
Thus, our final result for the spontaneous magnetization at the infinite volume limit is 



M* =0 = 0.710536 t 1 x 10" 6 (6) 
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-iog lo (/0- &) 


Gf(0) 


G|(0) 


G3 


1 


0.710537 


5.1449464 


-451.801295 


2 


0.372929 


147.744937 


-483052.394150 


3 


0.181173 


3270.020749 


-4.419043E8 


4 


0.0864639 


67534.150592 


-3.817470 x 10 11 


5 


0.0410479 


1.362846 x 10 6 


-3.231856 x 10 14 


6 


0.0194518 


2.727356 x 10 7 


-2.720537 x 10 17 


7 


0.00921183 


5.441076 x 10 8 


-2.280423 x 10 20 


8 


0.00436138 


1.084248 x 10 10 


-1.914112 x 10 23 


9 


0.00206473 


2.159631 x 10 11 


-1.597785 x 10 26 


10 


0.000977434 


4.301016 x 10 12 


-1.343532 x 10 29 


11 


0.000462716 


8.564077 x 10 13 


-1.113270 x 10 32 


12 


0.000219084 


1.704184 x 10 15 


-9.407884 x 10 34 



1 




UO 80-60 90,60 

jf uox^iezT^auSie]^ 



